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We derive the integrability conditions of nonautonomous nonlinear Schrodinger equations using 
the Lax Pair and Similarity Transformation methods. We present a comparative analysis of these 
integrability conditions with those of the Painleve method. We show that while the Painleve in- 
tegrability conditions restrict the dispersion, nonlinearity, and dissipation/gain coefficients to be 
space-independent and the external potential to be only a quadratic function of position, the Lax 
Pair and the Similarity Transformation methods allow for space-dependent coefficients and an ex- 
ternal potential that is not restricted to the quadratic form. The integrability conditions of the 
Painleve method are retrieved as a special case of our general integrability conditions. We also 
derive the integrability conditions of nonautonomous nonlinear Schrodinger equations for two- and 
three-spacial dimensions. 
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I. INTRODUCTION 



The question of integrability of nonautonomous nonlinear Schrodinger (NLS) equations has recently been extensively 
addressed due to their applications in the fields of trapped Bose-Einstein condensates and optical solitons in fibers 
Among many works that investigate the integrability of the celebrated NLS equation pJ-Q, He et al. [§] address 
the most general form of this equation, namely 

f(x, t) V m (x, t) + g(x, t) \9{x, t)\ 2 ^(x, t) + v{x, t) V(x, t) + i j(x, t) V(x, t) + i V t (x, t) = 0, (1) 

where the dispersion coefficient f(x, t), nonlinearity coefficient g(x, t), gain/loss coefficient 7(0;, t), and external poten- 
tial v(x,t) are real functions. Applying the Painleve test, the authors were able to derive the following integrability 
conditions 

f{x,t)=f{t), g(x,t)=g(t), j(x,t)= 7 (t), (2) 

v(x,t) =v (t) + v 1 (t)x + v 2 (t)x 2 , (3) 
where vq (t) and v\ (t) are arbitrary and vi (t) is given by 

4/V v 2 + .Dr . hi + fg) + g 2 (f 2 - //) - V J ir = 0. (4) 

Here and throughout, the dot represents time derivative and the subscripts represent partial derivatives. Clearly, 
this shows that, according to the Painleve test, Eq. ([T]) is integrable only for time-dependent dispersion, nonlinearity, 
gain/loss coefficients and for quadratic external potential. 

It is well-known, however, that an equation can be integrable in one sense and nonintegrable in another sense. 
Therefore, we use in this paper two different methods to investigate the integrability of Eq. ([1]). We use first the 
Lax Pair method and then the Similarity Transformation method to show that Eq. ([I]) is indeed integrable for more 
general time- and space-dependent coefficients. We derive an integrability condition that reproduces the conditions 
©-([I]) as a special case. It is found that, while the three methods generate different integrability conditions, when 
the coefficients of Eq. ([T]) are space- and time-dependent, the three methods generate the same integrability condition 
when the coefficients are only time-dependent. 

While it is not the aim of this paper to discuss and derive exact solutions of specific examples of integrable NLS 
equations, we present in section [IT] some examples corresponding to known integrable cases. In addition, we present in 
sections HT1 and Hill as a counter example to the Painleve result, two cases of NLS equations that are indeed integrable 
with space- and time-dependent coefficients. 

Finally, we use the Lax Pair method, in section HVl to generalize the integrability conditions for NLS equations of 
higher spacial dimensions. To the best of our knowledge, the integrability conditions derived here (Eqs. (J7J and (|23p) 
are presented in the literature for the first time. 
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II. LAX PAIR METHOD 



In the Lax Pair method for solving nonlinear partial differential equations, a pair of 2 x 2 matrices, U and V, is 
first constructed (or derived). The Lax Pair is used to define a linear system of equations for an auxiliary field <I>, 
namely <& x — U • $ and 3> t = V • <fr. The compatibility condition of this system, \J t — ~V X + [U,V] = 0, is required 
to be identical to the nonlinear differential equation. The Darboux transformation [lCf is then applied on the linear 
system to obtain new solutions from known solutions. 

We consider the NLS equation to be nonintegrable if the Lax Pair does not exist. Our method of Lax Pair search 
provides an answer to the question of existence of Lax Pair. In this method, we expand the Lax Pair U and V 
in powers of \& and its derivatives with unknown function coefficients and require the compatibility condition to be 
equivalent to the nonlinear differential equation. This results in a set of coupled differential equations for the unknown 
coefficients. Solving this system of equations determines the Lax Pair and results in an integrability condition on 
the functions f(x,t), g(x,t), j(x,t), and v(x,t). The details of this calculation are relegated to Appendix [SJ The 
integrability conditions of Eq. (fTJ) turn out to be 

7(M) = *M (6) 

g(x,t) 2c 2 [t) 

fg 3 (ft (g t - 257) - fug) + f?g 4 + 2/V (gv xx ~ g x v x ) 
+ f 2 g 2 (g (4<7*7 + gu) - 2 5t 2 - 2 g 2 ( 7t + 2 7 2 )) 

+ f 4 fag* - A8gg xx g 2 x + 10g 2 g xxx g x + g 2 [&g 2 xx - gg^)) = 0. (7) 

This condition shows that f(x,t), "/(x,t) and v(x,t) are determined by g(x,t) and the arbitrary functions C\(t) 
and C2(t). To get specific forms of integrable NLS equations, one starts with a certain form of g(x,t), Ci(t), and 
C2(t) from which the functions f(x,t), j(x,t), and v(x,t) will be determined according to Eqs. ([5][7]). With such 
a combination of coefficients, Eq. ((T|) will be integrable. Equations ([5][7|) represent a main result of the present 
paper. All previously-known special cases can be derived from these equations. We list here some of these special cases. 

Special case I: Constant and linear external potential 

With the choices: g(x,t) = 1 and c\(t) — C2(t) — 1, Eq. (TT]) takes the form 

+ |*| 2 * + f ca(i) + c 4 (t) x) * + i tf t = 0, (8) 



and 



where ca(t) and 04(f) are arbitrary real functions arising from integrating Eq. ([7]) with respect to x. Clearly, with 
= Ci(t) = 0, the well-known homogeneous Gross-Pitaevskii equation (GPE) is obtained. It is also established 
that Eq. © is integrable and exact solutions have been derived [l2[ • 

Special case II: Harmonic potential and gain/damping term 

For g(x,t) = 1, ci(i) = 1, and 02(f) = e at , where a is a real constant, Eq. ([I} takes the form 

^xx + |*| 2 * + x 2 - I i + c 3 (t) + ci(t) x^j * + % * t = 0. (9) 

With Ca(t) = c^f) — this will be the typical GPE with an expulsive harmonic potential and damping. 

Special case III: Harmonic potential and time-dependent nonlinearity 

For g(x,t) =e at , a(t) = c 2 (t) = e 2at , Eq. ||IJ) takes the form 

*xx + e at |*| 2 * + x 2 + c 3 (t) + c 4 (t) x^j * + i * t = 0. (10) 

With c 3 (t) = Ci(t) — this will be the GPE with harmonic potential and nonlinearity growing exponentially with 
time. This equation was shown to be integrable and exact solitonic solution were obtained by Ref. [J]. 
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Special case IV: ^"-dependent coefficients 

For g(x,t) = x n , where n is an integer, and c\(t) = c>}(t) = 1, Eq. (Q} takes the form 

x~ 2n V xx + x n |*| 2 * + (--n(n + 2)x- 2( - n+1 '> + x n+1 + cAt) x] * + i * t = 0. (11) 

\ 4 n + 1 / 

This is our first counter example to the conclusion of the Painleve test; an integrable NLS equation with space- 
dependent coefficients. 

Special case V: Time-dependent coefficients 

For the special case of g(x,t) = g(t), Eqs. (0 and dH) result in f(x,t) — f(t) and "f(x,t) — j(t), where f(t) and 7 ft) 
are now independent. In this case, Eq. (|7|) leads to 

-3 2 / 2 + fg (gf + f (2<?7 - <?)) + f 2 (2.g 2 -g(g + ^g) + 2g 2 (7 + 2 7 2 )) 
< x ^) = 4^2 x 

+ c 3 (t)x + Ci (t), (12) 



where 03(f) and c&(t) are arbitrary. Comparing this expression with v(x, t) — V2(t) x 2 + v%(t) x + V3, we get 

- Jtrl + fgfg - Vg 2 if + g 2 f 2 + fgg + Klin) - -'./-'//- 

+ 4v 2 fg 2 - 2/V 7 - 4/ 2 5 2 7 2 = 0. (13) 
This is the integrability condition (22) obtained by He et al. [§]. For the special case of 7(f) = 0, we get 

4/V v 2 + fg(fg + fg) + g 2 (f 2 - ff) - 2f 2 g 2 = 0, (14) 
which is the integrability condition (26) of He et al. 

III. SIMILARITY TRANSFORMATION METHOD 

This is a method that is used frequently in the literature to generate new solutions of a nonautonomous NLS 
equation from those of the standard homogeneous NLS equation 0-0j 03 • By transforming the coordinates and the 
wave function such that the nonautonomous NLS equation is transformed into the standard NLS equation, one can 
use the same transformation to obtain solutions of the nonautonomous NLS from the solutions of the standard NLS. 
Here, we exploit this method to investigate the integrability of NLS equations. It turns out that for the transformation 
to be possible, the coefficients of the nonautonomous NLS must satisfy certain integrability conditions. 

The integrability condition (4) can be obtained by requiring the transformation 

*0, t) = exp ( P(x, t) + i 0(x, t)) Q{X(x, t)), (15) 

where /3(x,t), 9(x,t), and X(x,t) being real functions, to transform Eq. ([TJ into the following time-independent 
homogeneous equation 

p(x,t) (eQ xx {X) + S\Q{X)\ 2 Q(X)^ -0, (16) 

where p(x,t) is in general a complex function that will be determined in terms of /3(x,t), 9(x,t), and X(x,t), and e 
and 5 are real constants. Substituting this form of ^(x, t) in Eq. ([T]), it takes the form 

e 20 g Q{X)\Q(X)\ 2 + f X 2 Q"(X) + [iX t + f (2X X {[3 X + i6 x ) + X xx )} Q'(X) 
v + i-i + ipt-Bt + f^ + ie^ + p^ + ie^ Q(X) = 0. (17) 



Comparing this equation with Eq. (|16|) . we get 8 coupled differential equations for the unknown functions. These 
equations can be solved resulting in rather lengthy expressions for the functions f(x,t), 9(x,t), and v(x,t) in terms 
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of /3(x, t) and arbitrary functions of t, as presented in detail in Appendix iBl We also show in Appendix [Bl that in the 
special case of /3(x,i) — f3(t) and j(x, t) — 0, the expression for v(x,t) simplifies to 



v(x,t) = c 5 



c 6 2 S 2 fc 2 



Ae 2 g 2 

ceS (fgc 2 + c 2 (gf - 2fg 

-fg (fg + fg) + <r (// - f 2 ) + V 2 g 2 
i/V 



,2 



(18) 



where cq is an arbitrary real constant, and C2{t) and c§{t) are arbitrary real functions. Equating the coefficient of the 
x 2 -term of the previous equation with V2, we retrieve the intcgrability condition Eq. (|4]) of the Painleve analysis and 
the Lax Pair method. 

As a second counter example on the integrability of Eq. ([1} with time- and space-dependent coefficients, we take 
the special case of (3(x, t) = —(1/2) log a: and 7(2;, t) = 0. Substituting this in f(x, t), g(x, t), and v(x, t), Eq. (TT]) takes 
the form 

i* i + * ax + t^3 W 2^ + Qg__y q , = Q (1Q) 

where we have set, for simplicity, 6 = e = 1. An exact solution of the homogeneous NLS, Eq. ([T5|. is 
Q(X) = — sn(X/v / 2 I — 1). Applying the Similarity Transformation, we get a solution of Eq. (Tl9|) . namely 
^/(x,t) = exp (—i x 2 /8 1) sn(t x 2 /V8 \ — l)/y/x, where sn(x\m) is the Jacobian elliptic function of modulus m. 



IV. INTEGRABILITY CONDITIONS OF TWO- AND THREE-DIMENSIONAL NLS EQUATIONS 

For the one-dimensional NLS equation, Eq. (fT]), to be integrable, the nonlinearity coefficient, g(x,t), and the 
dispersion coefficient, f(x,t), should be related according to the integrability condition Eq. ([5]). It is expected that 
an additional term proportional to ^ x { x i t) m t ne NLS equation would affect this relation between g(x, t) and /(a:, t), 
such that a wider class of integrable NLS would be obtained. In other words, it is hoped that one-dimensional NLS 
equations which are nonintegrable to become integrable as a result of this addition. In addition, with this term the 
NLS corresponds to two- and three-dimensional physical systems. Furthermore, such an addition would allow for the 
consideration of position-dependent effective mass situations. 

We choose the Lax Pair method to investigate the integrability of the following general NLS equation 

f(x, t) V xx (x, t) + h(x, t) ^ x (x, t) + g(x, t) |*(x, t)\ 2 ^(x, t) + v(x, t) *(x, t) + i j(x, t) &(x, t) + i ^ t (x, t) = 0, (20) 

where h(x,t) is a real function. 

Employing the Lax Pair method in a similar manner as explained in section [TT] and Appendix [XJ we obtain the 
following integrability conditions 

d(t)H(x,t) 

fM= 9 ( X ,tr ' (21) 

/ .x _ 9t{x,t) lc 2 (t) \ H t {x,t) 
1(,) ~g(x,t) 2c 2 (t) 4H(x,t)' 1 ' 

and 



±fg A H z (f t (2 57 - g t ) + f tt g) + if'g^H' + if g H (v x (gH x - 2g x H) + 2gHv xx ) 



+ f 



- g 2 H x (3 (gg xx - 2g 2 x ) H x + gg x H xx ) - 2gH (4&g*H x - lQgg 2 x H xx + gg x (gH xxx - A2g xx H x ) 

2g 2 (2g xx H xx + 3g xxx H x )) + 4 (365* - 48gg xx g 2 x + Wg 2 g xxx g x + g 2 (&g 2 xx ~ ,g 5 (4) )) H 2 
4f 2 g 2 H 2 [g (ig tl + g u ) - 2g 2 - 2g 2 ( 7t + 2 7 2 ))] = 0, (23) 
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where 

^ (M)=exp (/^4 (24) 

As a first check, we confirm that the above integrability conditions indeed reduce to those of the previous sections, 
Eqs. ([5][7]) for H(x,t) = 1 or equivalently h(x,t) = 0. These integrability conditions give all coefficients of the NLS 
equation in terms of two arbitrary functions g(x,t) and h(x,t). In the following, we present some interesting special 
cases. 

Special case I: Constant effective mass 

For g 2 (x,t) — H(x,t) = x n , we get the following integrable NLS equation 



i * t (a:, t) + ^(at, t) + - V x (x, t) + x n \*(x, t)\ J V{x, t) + I c 3 (t) + Ci(i) a: + ^ y j i) = 0, (25) 

where, 03(4) and C4(t) are arbitrary real functions. This equation corresponds to a NLS in n + I spacial dimensions 
with constant effective mass. Of particular importance, are the cases of n = 1 and n = 2, corresponding to two and 
three dimensions, respectively, with power law for the strength of the interatomic interaction, x n , and an effective 
potential that includes the centripetal part with a suitably-defined angular momentum. 

Special case II: Power-law effective mass 

For g 2 (x,t) = x p and H(x,t) = x q , the integrable NLS takes the form 

i f) + x^ 2 ^ xx {x, t) + \qx^- x * x (x, t) + x p \9(x, t)\H(x, t) 



(c 3 (t) + 2 l C *®_ q * 1+p - q/2 -\ P (2+p-q) I""*" 2 ) t) = 0. 



(26) 



The position-dependent effective mass appears in the NLS equation through the second derivative operator, which in 
this case will be given explicitly by 

d( 1 «(M)\_ i9 - 2 p $ji(M)+ l x«-*P-^ x {x,t). (27) 



dx \2m(x) dx J '2 
This can be satisfied for q = 4p, which corresponds to an effective mass m (x) = -x~ 2p /2, and Eq. (|55J| simplifies to 

iV t {x,t) +x 2pi & xx (x,t) +2 P x 2p - 1 ^ x (x,t) +x p \^(x,t)\ 2 V(x,t) 
+ U(t) + ^^x l - p -^p(2-3p)x 2 ^- 1 ^(x,t)=0. (28) 

V. CONCLUSIONS AND OUTLOOK 

Using the Lax Pair and Similarity Transformation methods, we have shown that, in contrast to the results of the 
Painleve analysis, Eq. ([T]) is integrable for time- and space-dependent coefficients. This is of course not surprising 
since it is well-known that an equation can be integrable in one sense and nonintegrable in another sense [l4j . For the 
special case when the coefficients of Eq. (JTJ are only time-dependent, the integrabilty conditions of the three methods 
become identical. 

The general integrability conditions found here, Eqs. and (|23|). would be useful for cases of space-dependent 
dispersion such as optical solitons propagating in a medium of space-dependent refractive index. 

It is in the nature of the Lax Pair method that new solutions are obtained using old (seed) solutions. Therefore, 
the integrability conditions imposed by the Lax Pair method are conditions on the possibility of mapping solutions 
of the nonautonomous NLS equation into other solutions of the same NLS. On the other hand, the integrability 
conditions imposed by the Similarity Transformation method are conditions on the possibility of mapping solutions 
of the homogeneous NLS equation into those of the nonautonomous NLS. The integrability conditions of the Painleve 
method have even more fundamentally different origin related to the movable singularities of the NLS equation. The 
fact that these three fundamentally different methods generate the same integrability condition for the case of time- 
dependent coefficients of the NLS equation, suggests two illuminating ideas which require further investigation. First, 
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the three seemingly fundamentally different methods may be after all not so different such that they can be unified 
in a single transformation method. Secondly, it seems that - at least for the case of time-dependent coefficients - 
integrability is an intrinsic property of the NLS equation, which is independent of the method used to reveal it. 

Finally, it is worth mentioning that all of the complicated calculations, including the Lax Pair search, simplify- 
ing, and solving the compatibility equations, that led to the results of this paper were facilitated using symbolic 
programming with the computer program Mathematica 7 [l5[ . 

Appendix A: Integrability conditions of the Lax Pair method 

The Lax Pair U and V are expanded in powers of fy(x,t) and its derivatives, as follows 

u _ ( «11 «12 ^ v _ ( «11 Wl2 

where 

Mil = fi(x,t)+f 2 (x,t)V(x,t), 
U12 = f 3 {x,t)+f 4 (x,t)V(x,t), 
«ai = f 5 (x,t) + f 6 (x,t)V*(x,t), 

U22 = f 7 (x,t)+fs(x,t)V*(x,t), 

v n =: 9i (x, t) +g 2 {x, t)H!(x,t) +g 3 (x,t)V x (x,t) + g 4 (x, t)V(x,t)V*(x,t), 
V12 =95(x,t) +g 6 (x,t)V(z,t) + g 7 (x, t)V x (x,t) + g 8 (x,t)^(x,t)^*(x,t), 
«ai = ga(x,t) + g 1Q (x,t)V*(x,t) + g n (x,t) %(x, t) + g 12 (x, t) V(x, t)V*(x, t), 
V22 =9i 3 (x, t) +g u (x,t)y*(x,t) +g 15 (x, t) V* x (x, t) + g 16 {x, t)V(x, t) t), 
and fi-s(x,t) and gi-%e{x,t) are unknown functions. 
The compatibility condition 

U,-V +\U VI -0- ( ° Pi(x,t)F[9] » 

U t V X + [U, Vj -U- y p ^ Xjt ) F *^ o )• (AD 

with pi (x, t) and p2{x,t) being unknown functions, requires 

h = 93 = fs = .915 = .98 = 3i2 = .914 = 32 = f 3 = fs = .99 = .95 = 0,/ 4 = ipi,f 6 = -ip2,gi = -fpi,9u 
-fP2,gi = -.9i6 = -ifViPi, 

and 



At-flix = 0> (A2) 

ht - 9n x = 0, (A3) 

2fpip 2 + 9 = 0, (A4) 

fxPi- fPi(h~ h) + fPi x ~ 96 = 0, (A5) 

/*P2 + /P2(A - fr) + fP2 x - .910 = 0, (A6) 

96 (fi - h) ~ ipi(gi ~ 913 - w + 7) - g 6x + ip u = 0, (A7) 
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Pio (/i ^ iV) + ipa(si - 5i3 - w - 7) + 510^ + ip2 t = 0, (A8) 



UpiP2)x + ffioPi +.96P2 = 0. (A9) 



1. Deriving a relation between f(x,t) and g(x,t) 



Substituting for fpiP2 from Eq. (IA4[) in Eq. (IA9|) . the latter simplifies to 



- ^.9x + 510 Pi + 96 P2 = 0. (A10) 
Multiplying Eq. (|A5[) by P2 and Eq. (|A6p by pi , and then adding, we get 

P2 (fPl)x+Pl(fP2)x -Pi 96 -Pi 910 
= {PlP2f)x+PlP2fx~P2 96 -Pl9lQ^0- (All) 

Substituting for P1P2, {f PiP2)xi and P2 96 +Piffio from Eqs. (|A4|) . (|A9|) and (|A10|) . respectively, the last equation 
takes the form 

k = -2^, (A12) 
/ 9 



which results in 



where c(t) is arbitrary. 



c(t) 

= ^55. (A13) 



2. Deriving a relation between w(a;, t) and g(x,t) 



Multiplying Eq. (|A5[) by P2 and Eq. (| A6|) by p\ , and then subtracting, we get 



Pi 

96 = — 
9 



h-h = -j^ + \Uo g ^ + l^. (A14 ) 
fpi 2dx g g 



(A15) 



Multiplying Eq. (|A7I) by p 2 arL d Eq. (|A8[) by pi, and then subtracting, we get 



£jfe ,iff (2c7 + c).g 
2ff 2 2 



ci + £ ku - + - I [3g t ; ) dx 



To obtain the last equation in this form, we have substituted for /1, gio, p%, and / from Eqs. (|A14|) . (|A10|) . (|A4|) . and 
(|A13[) . respectively. Here Ci(t), ki r (t), and ku(t) are arbitrary real functions that arise from two integrations over x. 
The last two being the real and imaginary parts of a complex function. The function c\(t) + k\ r {t) appearing in the 
last equation represents another arbitrary constant. Therefore, we can set, without loss of generality, k\ r (t) = 0. 
Multiplying Eq. (|A7[> by P2 and Eq. (|A8|) by p\ , and then adding, we get 



. t-9 96 , c , 2 9Pu + 2ig 6 g x - 3pi g t ic{gg xx -2g 2 x ) 
cp{ 2c 2gp x 2g 4 



To obtain the last equation in this form, we have substituted for f±, gio, P2, and / from Eqs. (|A14[) , (|A10[) . (|A4I) . and 
(|A13|) . respectively. Subtracting Eq. (|A3[) from Eq. (| A2|) . and substituting for /1 — /V and g\ — 513 from Eqs. (|A14|) 
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and (|A16I) . we get 



, c d 3 1 2k u g t g ■ 

VxM = ~2-gd^-g + — -c {2kul + ku) 

"' "~" n 1 (3gt - ~(2c 7 + c)) 



g 

2c 3 



(c (g t (2c7 + c) - 3cg M ) + ,g (c (2c7 t + c) - c 2 )) dx 

-(2017 + 0!)--^). (A17) 
c c / 

Since u is assumed to be real, the imaginary part of the last equation must vanish. This results in the following 
condition on 7 

(2ci 7 + ci) g-2c igt = 0, (A18) 

which gives 

^) = t|-^ (A19) 

Substituting this expression for 7(2;, t) in Eq. (|A17|) . results in following integrability condition for v x (x,t) 
v x (x,t) = n J 2 c 2 cici J - + 9tj dx 



2c 3 d 2 

(ci 2 (c(c5 t - C3 M ) +g(cc- c 2 )) + c 2 c\g - c 2 ci(cig + Ci^t)) dx 

c d 3 1 _ gjcjkn - kuCj) 
2 g dx 3 g cc\ 



(A20) 



This condition gives v x (x,t) in terms of c(t), ci(t) and g(x,t). Using Eqs. (|A13[) and (|A19I) to substitute for c(t) 
and ci(t), the last equation gives v x (x,i) in terms of f(x,t), 7(2:, t), g(x,t) and the arbitrary function knit) 

v x (x, t) = ^-j, 1 5 4 (/ ( 5t (/ t + 2/7) - / 5if ) + g (/ (/„ + 2/ 7t ) - / 2 )) 

+ 2/V (g t - 57) y («?(-| - 2 7 ) + fft ) 

1, ^ 1 2g(2fcn7 + fcn)~4fc Hgt 
" 2 /5 ^ ~g 27? ■ (A21) 

Combining the two integrals in Eq. (|A20|1 . differentiating with respect to x, and substituting for c and c\ in terms of 
/, g and 7, the last condition takes the form 

fg 3 (f t (g t - 2 57 ) - f tt g) + f t g A + 2/V (gv xx - g x v x ) 

+ fg 2 (g (4<?t7 + gu) - ig\ - %g 2 (n + 2 7 2 )) 

+ / 4 (36 5 4 - 48 5feff 2 + Wg 2 g xxx9x + g 2 (&g 2 xx - gg^)) = 0, (A22) 
which is identical to Eq. ([7]). 

To summarize: We have derived three integrability conditions: 

1. Condition (|A13|) gives f(x,t) in terms of g(x,i) and an arbitrary function c(t). 

2. Condition (|A19[) gives j(x,t) in terms of g(x,t) and an arbitrary function c\(t). 

3. Condition (|A22[) gives v x (x,t) in terms of g(x,t), f(x,t) and -f(x,t). 

The three arbitrary functions c(t), c\(t), and ku(t) are independent from each other. 
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Appendix B: Integrability conditions of the Similarity Transformation method 

Comparing Eq. (fT6j) with Eq. (fT7|) we get 

e 20 9 = 6p, (Bl) 

fX 2 x =ep, (B2) 

X t + 2fX x x = O, (B3) 

2X X f3 x + X xx = 0, (B4) 

/8t + 7 + /(2jM* + 0x*) = O, (B5) 



t + ^ + A«) - 0. (B6) 



Solving Eq. flB4j| for X, we find 



Eliminating p(x,t) from Eqs. (jBll) and (|B2[) and substituting for i) from the previous equation, we find a 
compatability condition that combines f(x,t) and g(x,i), namely 



X(x, t) = d(t) + C2(t) / e^ 2 * 3 ^ dar. (B7) 



ff (a;,t) = ^ C2 (t) 2 e - 6 ^)/(a ; ,t). (B8) 



Solving Eq. (|B3j for 0(a;,i), we get 



6(x, t) = - ^ or " dx + <*(*)> ( B9 ) 

J 2c 2 / 

where cg(i) is arbitrary. Solving Eq. (|B5I) for f(x,t), we get 

/ / ^ (4e^fe (J (eg - 2cgA) dx + cQ - 4c 2 /3 t + c 2 - 2 c 2 7 ) \ 

where c/(£) is arbitrary. 

1. Special case I 

To obtain the special case of time-dependent coefficients of the NLS equation, namely f(x,t) = f(t), g(x,t) = g[t) 
, we assume f3(x,t) = (3(t) and 7(x, i) = j(t). With these assumptions, Eq. (|B10|) takes the form 

2c 2 (2i, + -,)-c 2 

f(x,t)= ( K (2c 2 b-c 2 )x-e 2b c 1 ) 2 ^ . (Bll) 

The function f(x,t) becomes a;- independent either by setting the exponent in the previous equation to zero or by 
setting the coefficient of x to zero. We notice that the coefficient of x is identical with the denominator of the exponent. 
Therefore, we take the first choice. Setting the exponent to zero and solving for c 2 , we get 

c 2 (t) =c 8 e 4h+2 ^ dt . (B12) 

Substituting this expression for c 2 (i) in Eqs. (|B10[) and (|B9I) . we get f(x,t) = c/(i) = f(t) and 

x (xb + ci e' 2( - b+ -f t dt ) /c s ) + x 2 1 
9(x, t) = c e - * . (B13) 
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Finally, the function v(x, t) can be obtained from Eq. (|B6|) 

-2(H-/7«tt) (4/ (6 + 7)+/) -/Cl 



e 

u(a;, i) = - 



2cs/ 2 

/ (& + 7) + / (-6 + 4 7 6 + 26 2 - 7 + 2 7 2 ) 



2/ 



2 



x 2 



a2 -4(6+/ 7 dt) 

+ 4CV +6e (B14) 
taking the special case 7 = and substituting from Eq. (|B8[) for b = (1/2) log (e g/c\ 8 /), the last equation gives 

, ^ • C S 2 S 2 fcl 

V ( X >V = <% + 4£ 2 g2 

c 8 ^(/5Ci+C! (5/ -2/5)) 
27/^ x 

/<?(/<? • .///) • <r(// / 2 ) • 2/^ 
3^2 



4/ 3 ff 



x 2 (B15) 



Equating the coefficient of the cc 2 -term with V2, we retrieve the integrability condition Eq. (jj) of the Painleve analysis 
and the Lax Pair method. 
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